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Abstract. We compare the computation times and precisions of three expan- 
sions for the gravitational potential. The evaluation of the series by Selvaggi et 
al (2008 Class. Quantum Grav. 25 015013) is time-consuming to be applied to 
the gravitational constant (G) or STEP experiments. 

PACS numbers: 04.80.-y, 02.60.-x 

We compare the computational efficiency of three published expansions for the 
gravitational potential (GP) of a finite, hollow, thick-walled cylinder with constant 
mass density p and reflection symmetry with respect to the plane z = 0. The 
expansions are those developed by Lockerbie et al , Schlamminger et al and Selvaggi 
et al [DEI [3]. 

If V denotes the GP/(27rGp) and (2L,a,b) the length, inner and outer radii of 
the cylinder, Lockerbie et al formula, which is valid outside the sphere (S) of radius 
y/b 2 + L 2 and center at the origin, writes in terms of Legendre polynomials Pm as 
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The GP, Vscht^, z), by Schlamminger et al formula is expressed in terms of the 
GP along the axis of symmetry V(0, z). If we introduce f(x) = x + \Jb 2 + x 2 and 
g(x) = x\Jb 2 + x 2 , then V(0 7 z) & Vsch^, z) are evaluated by 

2V(0, z) = 6 2 {ln[/(z + L)] - ln[/(z - L))} + g(z + L) - g(z - L) - (a -> 6) , (2) 
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V Sch (r, z) = V(0, z) + J2 2 2 n(n!)2 ' ^ > (Converges for small r) . (3) 
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where (a — * b) substitutes a for b in the first four terms of (E|). 

In terms of the hypergeometric function 2-F1 @]j Selvaggi et al formula reads 

G(r, u, n, a, b) = b 2n+2 2 F 1 (n + 1, 2n + (1/2); n + 2; -b 2 /(u 2 + r 2 )j - (a -»■ b) , 
V ( (4n-l)ll fr\ 2n f z+L G(r,u,n,a,b) 
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We consider the case where a = 1/2, b = 1 & i = 1, in SI units. We restrict 
ourselves to the central plane z = and compare the three formulae (TTJ), & (HJ) 
in their respective regions of convergence. We choose three points: Pi (0.2,0) near 
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the z-axis, i-2(0-4, 0) near the inner surface of the cylinder and £3(2, 0) outside the 
cylinder. If V(P,n + 1) denotes the potential at the point P with n + 1 being the 
number of terms used (here, n represents the index of summation which has been 
given the same symbol in all formulae (P), & (JH)), our results follow where the 
exact decimal places are underlined 

V Se i{Pi,4) =i 0.820253 ; Vs e i(Pi, 7) ~ 0. 82038 2 ; Vs ch (Pi , 3) ~ 0. 820385 . (5) 
V Sci (P 2 , 5) ~ 0.825570 ; V Scl (P 2 , 7) ~ 0.828583 ; V Sch (P 2 , 5) ~ 0. 831536 . (6) 
V Scl (P 3 , 7) ~ 0.371846 ; V Loc (P 3 , 7) ~ 0. 37195 8 . (7) 

Vsch(-f3) blows up and so do Vl oc (Pi) & ^00(^2)- The number of exact decimal places 
in ([5]) to ([7]) has been identified upon comparing Vs c h(Pi, 3), Vs c h(P2, 5) & Vl oc (P 3 , 7) 
to more accurate results, which have been obtained at the same points using large 
values of n + 1 (see last paragraph). The relative errors at P are i?Sei-Sch(Pi, 4) — 
1.61 x 10- 4 , £ S ei-Sch(P2,7) ~ 3.55 x 10~ 3 & P S ei_Loc(P 3 , 7) ~ 3.01 x 10" 4 with 
£W-Sch(P, n + 1) = |F Se i(P, n + 1) - V Sc h(P, n + l)|/Vsch(P » + 1), etc. 

Near the plane z = 0, the expansions (p} & |3j) perform better than (j4|) (Eqs (HI 
to ([7])) and are easily handled with reasonable evaluation times. The evaluation 
of each term in ^ involves a non-trivial integration which is time-consuming. In 
contrast, dTJ) (J3j) are evaluated via the use of Legendre polynomials, which are built- 
in functions in most computer algebra systems, or by differentiation. For instance, the 
evaluation of Vs c i(P2, 7), by numerical integration, lasted more than 29000 times that 
of Vsch(P2,7) (if we re-derive ((2|) using computer-algebra, this ratio reduces to 370) 
and the evaluation of Vs e i(P3, 7) lasted more than 46000 times that of Vl oc (P3, 7). 

In the plane z = each term of ((4)) (n > 1 up to at least 5) reaches a maximum 
value for a < r^ < b and decays slowly near and off the axis resulting in a slow 
convergence of (j4]) around tq. Off the plane z = 0, the convergence of (j4]) improves 
around the axis due to an off-axis shift of ro and a drop in the maximum resulting 
in a rapid decay of its terms (n > 1) there. At the point P4(0.4, 2), we found 
V Sel (P 4) 7) = 0. 36482771 . Vs c h(Pi,7) = 0. 36482771 k V hoc (P 4 ,7) = 0. 36482 807. The 
evaluation of Vs c i(P4, 7) lasted 8500 times that of Vsch(Pi! 7) and 52000 times that of 

^Loc(P 4 ,7). 

The flexibility in using ((TJ or (J3j) allowed as to evaluate the GP using up to forty- 
one terms and to check over the stability of the results shown in ([5]) to (JTJ) for Vsch 
& VLoc- Moreover, the GP at each of the points P\ to P4 has been reobtained by 
numerically integrating the volume integral of the potential element di>'/(27r|r — r'|). 
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